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Abstract 

The reaction e"'"e^ e^e^7*7* e^e" hadrons for quasi-real photons is studied 
using data from a/s = 183 GeV up to 202 GeV. Results on the total cross sections 
(j(e+e~ — >■ e^e~ hadrons) and 17(77 ~^ hadrons) are given for the two-photon centre- 
of-mass energies 5 GeV < IV^^ < 185 GeV. The total cross section of two real 
photons is described by a Regge parametrisation. We observe a steeper rise with 
the two-photon centre-of-mass energy as compared to the hadron-hadron and the 
photon-proton cross sections. The data are also compared to the expectations of 
different theoretical models. 
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1 Introduction 



At high centre-of-mass energies, \fs, the two-photon process e+e~— >■ e+e^7*7*-^ e^e^hadrons 
is a copious source of hadron production. Most of the initial energy is taken by the scattered 
electron^. As their scattering angle is low, they often go undetected. The hadron system 
has, typically, a low mass compared to ^/s. A large fraction of the hadrons escape detection, 
due to the Lorentz boost of the 77 system and to the large diffractive cross section producing 
hadrons at small polar angles, where the detector acceptance is limited. For these reasons, the 
measured visible mass, W^is, is less than the two photon effective mass, W-y^. 

In this paper we analyse only data where the scattered electrons are not detected. New 
results on total cross sections (T(e"'"e~— > e'^e~ hadrons) are presented, using data collected with 
the L3 detector P for a total integrated luminosity of 51.4 pb~^ at \/s= 183 GeV, 171.8 pb~^ 
at ys= 189 GeV and 220.8 pb"^ at yi= 192, 196, 200, 202 GeV. 

The two-photon cross section cr (77 hadrons) is derived in the interval 5 GeV < W^^ < 
185 GeV, while the analysis of the data taken at ■\/s= 133 GeV and 161 GeV covered only 
the interval 5 GeV < W^^ < 75 GeV. 

2 Measurement of cross sections 

2.1 Monte Carlo simulation 

The e"'"e~— e"'"e~7*7*— > e~^e~ hadrons processes are generated with the PHOJET |2] and PYTHIA 
ji] event generators. For the annihilation process e"'"e~— > qq(7), PYTHIA is used. KORALZ jH] 
is used for e"'"e^— r+r~(7) and KORALW B] for e+e"— W"^W~. The e"'"e~— e'^e~T~^T~ chan- 
nel is generated by DIAG36 "7 . The Monte Carlo events are simulated in the L3 detector using 
the GEANT |3 and GHEISHA jU] programs and passed through the same reconstruction pro- 
gram as the data. Time dependent detector inefficiencies, as monitored during the data taking 
period, are also simulated. 

2.2 Event selection 

The analysis is based on the central tracking system, the high resolution electromagnetic 
calorimeter, the hadron calorimeter and the luminosity monitor. 

Two-photon events are collected predominantly by the track triggers ^Ij. The trigger 
efficiency is studied separately for each data sample by comparing the number of events accepted 
by the track trigger and the calorimetric energy trigger. The efficiencies of higher level triggers 
are measured using prescaled events. The trigger efficiency increases from 80 % at Wvis = 5 GeV 
to 94 % above 80 GeV. 

Hadronic two-photon events are selected by the following criteria : 

• To exclude scattered electrons, events with clusters in the luminosity monitor having en- 
ergy greater than 30 GeV, in a fiducial region of 33 mrad < < 64 mrad are rejected. The 
virtuality of the interacting photons, Q^, is thus less than 8 GeV^, with an average value 
< >^ 1.5 X 10-2 GeVl The distribution of low energy clusters in the luminosity mon- 
itor, presented in Figure shows a good agreement with both Monte Carlo programs. 
When the scattered electron reaches the detector, the agreement is maintained with the 

^■'Electron stands for electron and positron throughout this paper. 
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PHOJET Monte Carlo, while these configurations are missing in PYTHIA because of a 
p— mass cutoff, < m^, applied to the two-photon luminosity function in the generation 
of the events. 

• The total energy in the electromagnetic calorimeter is required to be greater than 0.5 
GeV, in order to suppress beam-gas and beam-wall backgrounds, and less than 50 GeV, to 
exclude radiative events, e+e~ — * Z7. The total energy deposited in the electromagnetic 
and hadronic calorimeters, E^ai, must be less than 40% of y/s, to exclude annihilation 
events, as shown in Figure^. 

• At least six particles must be detected, in order to exclude events containing r. A particle 
is defined as either a track, a photon in the electromagnetic calorimeter, or a cluster 
in the hadron calorimeter or in the luminosity monitor. Clusters in the luminosity mon- 
itor are considered as pions if their energy is below 5 GeV, as photons otherwise. The 
distribution of the number of particles is presented in Figure This distribution is not 
well reproduced by the Monte Carlo simulations. 

After selection, the background from beam-gas and beam-wall interactions is found to be neg- 
ligible. The visible effective mass of the event, VFvis, is calculated from the four-momenta of 
the measured particles. The analysis is limited to events with W^^g > 5 GeV. 

Almost 2 million events are selected, 1.6 x 10^ at a/s= 183 GeV, 7.8 x 10^ at y/s= 189 GeV 
and 1 X 10^ at \/s= 192 — 202 GeV. The average centre-of-mass energy of this last sample is 
v/i= 198 GeV. 

The W^vis spectrum is shown in Figure ^ for the total data sample. The background is 
below 1% at low masses, where it is dominated by two-photon r-pair production. It increases 
at high masses, due mainly to annihilation processes and reaches a maximum of 15%. 

The distributions of the rapidity, y, of the particles and of their energy flow are compared 
to the Monte Carlo expectations in Figure El A good agreement is observed also in the regions 
where | ?/ | ~ 3, between the luminosity monitor and the hadron calorimeter. 

2.3 Unfolding and efficiency 

The distribution of the two-photon effective mass W-y^ is obtained from the visible effective 
mass Wvis by the same unfolding procedure ^] used in Reference El For each data sample, 
the W^is spectrum is subdivided in 16 intervals, presented in Figure and the resulting 
Wryy distribution in 8 intervals, presented in Figure Eb- The result of the unfolding procedure 
depends on the Monte Carlo used. Data unfolded with PYTHIA are in general higher than 
if unfolded with PHOJET. After unfolding, the events are corrected for the efficiency, using 
the ratio between selected and generated events in each interval. This includes the purely 
geometrical acceptance as well as the efficiencies of the detector and the analysis procedure. 
For W^^> 30 GeV, the efficiency is rather constant, with a value of about 80%. The efficiency 
obtained with PYTHIA is lower by about 10%, which may be attributed to a different modeling 
of the diffractive interactions, of difficult detection. 

2.4 Cross Section Determinations 

The measured cross sections Acr(e+e^^ e~^e~ hadrons) are given in Table 1 for the three data 
sets, as a function of the W^^ intervals. The average of the results obtained by unfolding the 
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data with PYTHIA or PHOJET is used. Due to the unfolding procedure, the measurements 
are highly correlated. The correlation matrix, similar for the three data sets, is given in Table 2 
for 189 GeV. The differential cross section Aa/AW^-y is shown in Figure EJ together with 
our measurements at lower LEP collision energies ,2j. The fast decrease of the cross section as 
a function of W-y^ is due to the two photon luminosity function, which depends on W^^/s. 

The systematic uncertainties are evaluated for each W^^ bin. They are independent of the 
data sample, inside statistical accuracy. They are evaluated as follows and their contribution 
is listed in Table 3. 

• Trigger efficiencies: by varying this quantity within the accuracy of its determination, of 
about 10%. 

• Energy scale of the calorimeters and contribution of the annihilation background: by 
varying the Eca\ cut by ± 10% of y/s. 

• Uncertainties on the rejection of scattered electrons: by changing the Eiumi cut from 30 
GeV to 50 GeV. 

• Uncertainties on the particle multiplicity: by accepting a minimum number of four or 
eight particles instead of six. 

• Uncertainties due to Monte Carlo statistics are negligible for W-y^ < 65 GeV, but impor- 
tant in the higher W^^ bins. 

Uncertainties on the energy scale of the small angle calorimeter, evaluated by varying the gain 
by a factor two, are negligible. The total experimental systematic uncertainty, obtained by 
adding in quadrature all contributions, is also given in Table 3. The uncertainty related to 
the Monte Carlo model is given in the last column of Table 3. It is half of the difference 
between the results obtained by unfolding the data with PHOJET or PYTHIA and exceeds 
the experimental uncertainty in almost all bins. 

To extract the total cross section of two real photons, the luminosity function C^^ jT^I is 
calculated and the hadronic two-photon process is extrapolated to = 0. This is done as in 
Reference j2] by using an analytical program ^Hj- Depending on the choice of photon form 
factors, this calculation varies of ±5%. 

The cross sections obtained for the three data sets are compatible within the experimental 
uncertainties and are presented in Table 4. As expected from the study of the experimental 
systematics, the largest differences are observed in the first and last bins. The combined value 
is also given in Table 4 and in Figure EK with the statistical uncertainties obtained from the 
unfolding and the experimental systematics. The values obtained by unfolding the data with 
the two Monte Carlo programs separately are shown in Figure Eb and can be obtained from 
the last column of Table 4. 

3 Comparison with Theoretical Models 
3.1 Regge parametrisation 

The total cross sections for hadron-hadron, cipp, photon-hadron, a^p, and photon-photon, a^-y, 
production of hadrons show a characteristic steep decrease in the region of low centre-of-mass 
energy, followed by a slow rise at high energies. From Regge theory pBJ this behaviour is 
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understood as the consequence of the exchange of Regge trajectories, a{t), in the t-channeL 
The total cross section takes the form (Xtot oc The low energy region is sensitive to the 

exchange of a Reggeon R [R = p, u, f, a ..), with ai?,(0) ~ 0.5. At high energies, the Pomeron 
exchange dominates, with ap(0) ~ 1. A parametrisation of the form 

a,,, = As' + Bs-^ (1) 

accounts for the energy behaviour of all hadronic and photoproduction total cross sections, the 
powers of s being universal [16j. This is confirmed by the recent compilation of the total cross 
section data [Hj where a fit of Equation 1 for all hadron total cross sections gives a result 
compatible with the universal values e = 0.093 ± 0.002 and rj = 0.358 ± 0.015. The coefficients 
A and B are process and dependent. If photons behave predominantly like hadrons, this 
expression may also be valid for the two-photon total hadronic cross section, with s = W^^. 

Considering only the experimental uncertainties, statistical and systematic, several Regge 
fits are performed on the data and their results are presented in Table 5. The exponent r] 
is fixed to the universal value, since the low mass range is too small to be sensitive to this 
parameter. When the W-^-^ interval is restricted to 5 GeV — 65 GeV, a range similar to the 
one covered by our previous data |2], similar values of the parameters A and B are obtained. 
In this limited interval the data are compatible with the universal value of e. Extending the 
range to the whole W^^ interval, the fit with the exponents e and r] fixed to the universal value, 
dashed line in Figure EK, does not represent the a^^ energy dependence. A fit with A, B and 
e as free parameters, represented as a full line in Figure EK, gives e = 0.225 ± 0.021 with a 
confidence level of 4%. This value is more than a factor two higher than the universal value. 
It is independent of the Monte Carlo model used to correct the data, as shown in Table 5 and 
in Figure Eb- 

The fitted value of e is strongly correlated to the Reggeon component. To avoid this corre- 
lation, we fit only the Pomeron exchange for sufficiently high W^^ values. The results, using a 
different initial value of IV-y^, are listed in the second part of Table 5 . The value of e increases 
by increasing the lower mass cutoff, thus indicating that its value is not universal, but it reflects 
the onset of QCD phenomena, as e increases with increasing PV^^. 

3.2 Models for 77 total cross sections 

Several models P^IEU] were recently compared to the L3 and OPAL |2I] measurements. Their 
predictions for the two-photon total cross section are typically derived from measurements of 
proton-proton and photoproduction total cross sections via the factorization relation: o"^^ ~ 
cr^p/cTpp P2]- In general, these models give an energy dependence of the cross section similar to 
the universal flt discussed above. Two examples ^011201 are shown in Figure in comparison 
with the results of previous experiments 23j, those presented in this letter and those of OPAL. 
While the measurements at the low energy colliders present a wide spread, a good agreement 
is found between L3 and OPAL in the common W^^ range, 10 GeV < W^^ < 110 GeV. Good 
agreement is also found if the data, unfolded separately with either PHOJET or PYTHIA, are 
compared. In this PV^-y region, a model ^2Qj reproduces well the data and the predictions of the 
other ||l9j are too high by 20%. However, for both lower and higher values of W^^, the L3 data 
show a much steeper energy dependence than the theoretical predictions. 

In the Regge theory, the Pomeron intercept is 1, yielding a constant total cross section. 
When the rise of the proton-proton total cross section was first observed, it was explained [21] 
with an increase of the number of hard partonic interactions. The predictions of a model [2^] 
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that calculates such effects, using an eikonalized prescription to enforce unitarity, are shown in 
Figure IHb. The parameters of the model are determined from photoproduction data and the 
L3 results are well inside the uncertainty related to this extrapolation. 

Models with two Pomerons were recently proposed [2^1 to explain the fast energy increase 
of charm production at HERA. In this model, the 'soft' and the 'hard' Pomeron have different 
intercepts. Because of the qq component in the photon wave-function, the 'hard' Pomeron 
can contribute to the two-photon cross section even at = 0. Thus a more rapid energy 
dependence for cr^^ is expected. The increase in e with larger values of W^^, as listed in the 
second part of Table 5, is consistent with such a contribution of the 'hard' Pomeron. 
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183 GeV 


189 GeV 


198 GeV 


A W^^(GeV) 


< W^-^ >(GeV) 


AcTe+e- (nb) 


AcTe+e- (nb) 


AcTe+e- (nb) 


5 


-9 


6.7 


5.145 ± 0.025 


4.996 ± 0.009 


5.093 ± 0.008 


9 


- 17 


12.3 


3.358 ± 0.013 


3.350 ± 0.006 


3.466 ± 0.005 


17 


-31 


22.7 


1.812 ± 0.007 


1.880 ± 0.004 


1.962 ± 0.003 


31 


-47 


37.8 


0.776 ± 0.004 


0.813 ± 0.002 


0.857 ± 0.002 


47 


-65 


54.8 


0.388 ± 0.003 


0.422 ± 0.002 


0.453 ± 0.002 


65 


- 105 


80.2 


0.308 ± 0.003 


0.353 ± 0.002 


0.386 ± 0.001 


105 


- 145 


120.4 


0.070 ± 0.001 


0.096 ± 0.001 


0.111 ± 0.001 


145 


- 185 


158.7 




0.021 ± 0.001 


0.028 ± 0.001 



Table 1: The measured cross sections A(T(e''"e~— > e~ hadrons) as a function of the 77 centre- 
of-mass energy W^^ for the three data sets. Only the statistical uncertainties, obtained after 
unfolding, are given. 



A W-,„,(GcV) 


5-9 


9-17 


17-31 


31-47 


47-65 


65-105 


105-145 


145-185 


5-9 


1. 
















9-17 


0.931 


1. 














17-31 


0.815 


0.939 


1. 












31-47 


0.692 


0.803 


0.908 


1. 










47-65 


0.525 


0.602 


0.689 


0.761 


1. 








65 - 105 


0.336 


0.384 


0.436 


0.497 


0.486 


1. 






105 - 145 


0.130 


0.150 


0.166 


0.186 


0.190 


0.208 


1. 




145 - 185 


0.063 


0.072 


0.077 


0.080 


0.077 


0.089 


0.094 


1. 



Table 2: The correlation matrix after unfolding, for the data taken at \/s— 189 GeV. 



A Vr^^(GeV) 


Trigger 


-E'cal 


-^lumi 


-^part 


MC Stat. 


Total exp. 


MG model 


5-9 


0.9 


0.1 


< 0.1 


7.4 


< 0.1 


7.5 


7.0 


9-17 


0.9 


0.1 


< 0.1 


5.0 


< 0.1 


5.1 


1.2 


17-31 


0.7 


0.1 


< 0.1 


3.2 


< 0.1 


3.3 


4.0 


31-47 


0.6 


0.2 


0.2 


2.0 


< 0.1 


2.1 


6.8 


47-65 


0.5 


0.2 


0.1 


1.7 


< 0.1 


1.8 


8.9 


65 - 105 


0.4 


0.2 


1.1 


1.7 


1.5 


2.6 


10.4 


105 - 145 


0.4 


0.3 


3.0 


1.3 


8.1 


8.7 


15.5 


145 - 185 


0.4 


0.8 


6.4 


2.2 


12.4 


14.1 


27.4 



Table 3: Evaluation of systematic uncertainties due to the trigger, the analysis cuts and the 
Monte Carlo statistics. All values are per-cent uncertainties on the cross sections A(7(e+e~— > 
e~ hadions) and (7(77 — > hadrons). The uncertainty introduced by unfolding the data with 
PYTHIA or PHOJET is considered separately in the last column. A further scale uncertainty 
of 5% must be added for the a{jj hadrons) cross sections, due to the two-photon luminosity 
function. 
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183 GeV 


189 GeV 


198 GeV 


all data 






<W^^> { GeV) 


77 V / 


77 V / 


77 V / 


a^^i nh) 

77 V / 


AaZ^ (nb) 


A(T^^' (nb) 

'7'7 V / 


6.7 


422.6 ± 4.0 


394.9 ± 0.7 


398.4 ± 0.6 


397.2 ± 0.5 


30 


=F 28 


12.3 


378.4 ± 2.8 


360.2 ± 0.7 


368.2 ± 0.6 


365.2 ± 0.4 


19 


=F 4 


22.7 


359.8 ± 2.9 


348.9 ± 0.7 


358.2 ± 0.6 


354.4 ± 0.5 


12 


± 14 


37.8 


382.1 ± 4.5 


368.4 ± 1.1 


379.2 ± 0.9 


374.8 ± 0.8 


8 


± 26 


54.8 


408.6 ± 6.4 


403.0 ± 1.8 


418.0 ± 1.5 


411.5 ± 1.1 


7 


± 37 


80.2 


461.2 ± 8.6 


459.5 ± 2.2 


478.5 ± 1.9 


470.3 ± 1.4 


13 


± 49 


120.4 


496 ± 19 


556.7 ± 5.3 


586.1 ± 4.3 


572.0 ± 3.3 


53 


± 89 


158.7 




726 ± 15 


738 ± 11 


734.1 ± 8.7 


102 


± 202 



Table 4: The o"(77 hadrons) cross sections as a function of the average 77 centre-of-mass en- 
ergy, < W^^ >, for the three data sets and for their combination. The statistical uncertainties, 
obtained after unfolding, are given for each data set. The experimental systematic uncertainty, 
A(t!^^p, and the difference, Aa^, between the average value and the result unfolded with PHO- 
JET (lower sign) and with PYTHIA (upper sign) are also given. A further scale uncertainty of 
5% must be added, due to the two-photon luminosity function. 



Unfolding 


W^^ interval 


A 


B 


T] fixed 


e 


xVd.o.f. 


C.L. 




5-65 GeV 


178 ± 5 


453 ± 101 


0.358 


0.093 fixed 


5.3/3 


0.15 




5 - 185 GeV 


181 ± 3 


321 ± 120 


0.358 


0.093 fixed 


55/6 


10-9 




5 - 185 GeV 


58 ± 10 


1020 ± 146 


0.358 


0.225 ± 0.021 


12/5 


0.04 


PHOJET 


5 - 185 GeV 


52 ± 11 


1201 ± 146 


0.358 


0.221 ± 0.023 


8.6/5 


0.12 


PYTHIA 


5 - 185 GeV 


63 ± 10 


842 ± 146 


0.358 


0.228 ± 0.018 


19/5 


0.002 




17-105 GeV 


165 ± 21 






0.116 ± 0.016 


4.3/2 


0.12 




31-185 GeV 


113 ± 19 






0.163 ± 0.021 


3.4/3 


0.33 




47 - 185 GeV 


81 ± 23 






0.202 ± 0.035 


1.5/2 


0.48 



Table 5: Fits to the total cross sections listed in Table 4 of the form a^^= As^ + B s'"^ [16j, 
where s = W^^ . PHOJET and PYTHIA indicates that only this Monte Carlo is used to unfold 
the data. In all other cases the average unfolding result of the two generators is used. The 
statistical and experimental uncertainties and the correlation matrix between the data points 
are used. The fitted parameters are strongly correlated. The second set of fits evaluates only 
the increase of cr^^ with s, i.e. the Pomeron part of the fit. The values of the ^ind the 
corresponding confidence level are given. 
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Figure 1: Example of the distributions used for the event selection: a) energy in the luminosity 
monitor; b) energy Ecai in the electromagnetic and hadronic calorimeters, normalised to ^/s^, c) 
number of particles measured in the detector, d) The distribution of the visible mass Wyis, for 
the full data sample at ^/s — 183 GeV — 202 GeV. The data are compared with Monte Carlo 
predictions. The backgrounds due to e+e~ annihilation and e+e"— > e+e~T+T~ are indicated as 
a shaded area. 
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Figure 2: a) The distribution and b) the average energy of the measured particles as a function 
of the rapidity y, for the full data sample at a/s = 183 GeV — 202 GeV. The data are compared 
with Monte Carlo predictions. The backgrounds due to e+e~ annihilation and e+e~— > e+e-r+r" 
are indicated as a shaded area. 
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Figure 3: Example of the unfolding for data at y/s = 189 GeV. a) Distribution of the measured 
visible mass W^is in the 16 bins used for the unfolding, b) The measured W^is and the resulting 
W^^ spectra, obtained by unfolding from W^is to W^^ with the Monte Carlo PHOJET or 
PYTHIA. 
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Figure 4: The cross section Aa{e^e~ — > e'^e~ hadrons) / AW^j measured at ^/s — 133 GeV — 
202 GeV. Statistical and systematic uncertainties are added in quadrature and are often smaller 
than the symbol size. 
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Figure 5: The two-photon total cross section from the combined data at ^/s = 183 GeV — 
202 GeV. a) The average result, obtained by unfolding the data with the two Monte Carlo 
models, is used. Two Regge fits, described in the text, are superimposed to the data. The 
continous line corresponds to the fit with the coefficient e left as a free parameter, the dashed 
line is the fit with e fixed to 0.093. b) The two-photon total cross section obtained by correcting 
the data sample with PHOJET (full points) and with PYTHIA (open points). The Regge fits of 
Table 5 are superimposed to the data. The statistical and experimental systematic uncertainties 
are added in quadrature. 
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Figure 6: The two-photon total cross sections compared to various models, a) The predictions 
of References Cni and 1201 are compared to all two-photon total cross section data j2Tl|2H] . b) 
Predictions of the minijet model |2^; the two lines correspond to different choices of the model 
parameters. The statistical and experimental systematic uncertainties are added in quadrature. 
The uncertainties due to Monte Carlo models and to the two-photon luminosity function are 
included in the dashed hues in b). 
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